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ANALYSIS OF THE DYNAMICS OF A NUTATING BODY 
by William J. Anderson 
Lewis Research Center 


SUMMARY 

The equations for the displacement, velocity, and acceleration of a point in a nu- 
tating body are developed . These are used to derive equations for the inertial moment 
developed in a nutating body of arbitrary shape. The moment equations were applied to 
calculate the inertial moment acting on the support bearing system in a nutating plate 
transmission. For calculation purposes a previously designed nutating plate transmis- 
sion was used. The inertial moment was found to be 28 times as great as the moment 
resulting from reaction forces. This indicates that nutating plate drives are severely 
speed limited. 

The equations developed are considered to be potentially useful for calculating 
squeeze film and distortion effects in thrust bearings and face seals which are subject to 
high-frequency, low -amplitude nutation. 


INTRODUCTION 

The development of mechanical power transmissions that utilize a nutating member 
(ref, 1) has led to a need for a thorough understandir^ of the dynamics of a nutating body. 
Figure 1 (from ref. 2) shows a sketch of the arrangement of the components in a split nu- 
tating drive main rotor gearbox. The nutating plates describe a nutation cycle for each 
revolution of the input shaft. The nature of the nutating motion considered here can be 
visualized by referring to figure 2(a). The X,Y,Z coordinate system is fixed in the nu- 
tating body. The X axis describes a cone of angle 2a in space. Points along the Y 
and Z axes describe lemniscate-like figures. There is no rotation of the body about the 
X axis. In this respect the nutation considered here differs from that of spinning bodies 
that rotate about the nutation axis. Pins on the nutating plates (fig. 1) engage teeth on the 
fixed reaction plates and also on the drive plates. The drive plates advance one tooth for 
each revolution of the input shaft. 



The loads on the nutating plate support bearing, which result from both contact and 
inertial forces, must be known. In reference 2 a design analysis showed that contact 
forces on the nutating plate could result in severe nutator support bearing loads and in 
short bearing life. An approximate analysis of inertia loads assuming simple harmonic 
motion of the nutating body revealed that these, too, could be very significant, especially 
at high speeds, 

A search of the literature revealed a number of papers dealing with the dynamics of 
nutating bodies such as references 3 to 6. All of these, however, are concerned with the 
force systems that create the nutation and damping techniques for attenuating the nutation. 
None of them treat the dynamic forces generated by the nutating body. 

Therefore, the object of the work reported herein was to develop the displacement, 
velocity, and acceleration equations for a point in a nutating body and to apply these to 
the calculation of inertial moments. The resulting equations can also be used in the de- 
sign of nutating transmissions since they describe the kinematics of nutating motion. 

They may also be useful in the design and performance prediction of seals since nutating 
motions may occur in face seals. High-frequency nutating motions, even of very small 
amplitude, could conceivably play an important role in the dynamics of seals. 
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SYMBOLS 

inertial moment, N* sec^/cm‘* (lb - secVin. 
radius, cm (in.) 
thickness, cm (in.) 
time, sec 

volume, cm*^ (in. ) 

coordinate system fixed in the nutating body 
fixed inertial coordinate system 
point coordinates in X,Y, Z coordinate system 
point coordinates in Xq,Yq,Zq coordinate system 
nutation angle (between X and Xq, fig. 2), deg 

projection of cp into the X^Y^ plane (fig, 2), arcsin(sin or cos wt/cos xf/), 
deg 

harmonic oscillation angle, deg 

projection of r}/ into the XqZq plane (f^. 2), arcsin(sin or sin u;t/cos (p) 
angular coordinate for integration purposes 
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p density, N* secVc®^ (lb- secVin.'^) 

<p angle in XYqY plane between X and intersection of XYqY and XqZq planes 
(fig. 2), arcsin(sin of cos ct?t) 

\l/ angle in XZZq plane between X and intersection of XZZq and XqYq planes 
(fig. 2), arcsin(sin a sin wt) 

w nutation frequency, rad/sec 
Subscripts: 
h harmonic 
I inertial 
i inner 
o outer 


ANALYSIS 

Consider a body of arbitrary shape subjected to a nutating motion at a constant fre- 
quency Cl} and nutation angle a. Displacements, velocities, and accelerations of a 
point P(xj,yj,Zj^) will be derived relative to a fixed set of inertial coordinates Xq,Yq, 
Zq. The nutation motion will be described by the body fixed coordinate system X, Y, Z , 
as shown in figure 2. The axis of nutation is the X^ axis. 

The nature of nutating motion is such that the Yq axis is always in the XY plane, 
and the Zq axis is always in the XZ plane. 


Displacement Equations 

Xq = cos Of - yj^ sin <p cos e + Zj cos /3 sin (la) 

yQ = yj cos + Xj sin a cos wt + Z|^ sin i|/ sin j3 (lb) 

Zq = Zj cos - Xj sin or sin o;t + y^ sin <p sin e (Ic) 

In terms of a and cot equations (1) become 
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Xq = Xj cos a 


Yj sin 2a cos cot 


2(1 - singer cos^cot) 2(1 - sin^a sin^ 


Zj sin 2of sin cot 


a sin cot) 


2 2 sin Of sin 2cot 

Yq = yi(l - sin a cos cot) + Xj sin a cos cot + = (2b) 

n n 1/2 

2(1 - sin'^a sin^cot) 


.„2 . 


2 2 ^1 2cot 

Zq = Zj(l - sin a sin^cot) - Xj sin a sin cot + 


2 2 1/2 
2(1 - sin^or cos"^cot) 


VelocitY Equations 

^ _ CO sin 2o! sin cot „ , co sin 2 q; cos cot /„ s 

*0 T'yi (3a) 

, .2 2^/2^ „ 9 3/2 ^ 

2(1 - sin a cos cut) 2(1 - sin'^of sin'^cut) 

y = _(co sin or sin cot)x, + co sin^a sin 2cot ^ ^ co sm^a(cos 2cot + sin^a sin^cot) ^ 

. 2 2 .1/2 1 2 9 3/2 1 

2(1 - sin'^o; cos'^cot) (1 - sin^o; sin'^cot) 


L = -(w sin a cos cot)x, + sin^o>(cos 2cot - sin^a cos^cot) _ co sin^a sin 2cot 

1 _ Q /9 1 1 /o 1 




2 9 3/2 

(1 - sin O' cos cot) 


2 9 1/2 

2(1 - sin^a sin'^cot) 


Acceleration Equations 


, ^sin 2a(cos - 3 sln^n> co s ^ 2 sin^a cos^c^t) . sin 2a(-sin ^ + 3 sin^a sin a - 2 .sin^,> 


2(1 - sin^t:^ cos^wt) 


p 9 5/2 

2(1 ^ sin a sin (ijt) 
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yn = -(o>2 sin a cos -sin^o-cos"^ 

22^/2^ 

(1 - sin « cos a>t) 

22 2^ 2 45\ 

^ u) sin^ai(-2 sin 2a;t + 2 sin O! sin wt cos wt + 3 sin a sin at cos at - sin a sin at cos at) 
+ — 

2 2^/2 
(1 - sin"^Qf sin"^at) 




4(b) 


2 9 2 ^ 2 4 ^ 

z = sin a sin wt)x + ^ sin ty(-2 sin 2ciJt + 2 sin"^a cos'^o^t sin ojt + 3 sln"^g cos cot sin g*t - sin a cos cjt sin at) 

c /o 1 


2 2 

(1 - sin a cos cot) 


4(c) 


2 2 2 4 s 

to sin cv(cos 2cot 4- sin sin co^t) , 


(1 - sin^» sin^ort;) 


3/2 


Inertial Moments 

The convention used to calculate inertial moments developed by an arbitrary nutating 
body is shown in figure 3. The moment equations then become 

^ Vo - 

^yo = £ ^VO - VO^^ 

= f ^Vo - Vo^p 

u 

RESULTS AND DISCUSSION 

For purposes of calculation, conditions approximating those of the transmission de- 
sign developed in reference 2 were used. This design (ref. 2) was for an input power 
1. 86 megawatts (2500 hp), and a 28:1 speed reduction. The nutating plate pitch diameter 
was 91. 5 centimeters (36 in. ). A point was chosen on the pitch diameter of the nu- 
tating plate (see fig. 4). Values of w = 650 radians per second and of = 9*^ were used 
for calculating displacements, velocities, and accelerations. 
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Displacements, Velocities, and Accelerations 


Figure 4 shows the generalized motion of Pj projected into both the X^-Zq and 
Xq-Yq planes. The displacements projected into the X^-Zq plane Is a lemniscate-like 
curve (fig. 4(a)). In the Xq-Yq plane the displacement describes an arc (fig. 4(b)). 
Locations of Pj at various values of u>t are shown in figures 4(a) and (b). In the case 
of the nutating plate transmission, these displacement curves are necessary for the 
proper design of the engagement between the nutator and the drive and reaction plates. 

Figures 4(c) and (d) show the velocities of point Pj^ projected into the Xq-Zq and 
XqYq planes, respectively. The velocities x^, y^, and Zq can be added vector ially to 
get the absolute velocity of Pj at any point in the nutation cycle. This is necessary for 
the calculation of the relative velocity of engagement between the nutator and the drive 
and reaction plates. This, in turn, is required for the calculation of elastohydrodynaraic 
film thicknesses, traction forces, temperat;ire rise, power loss, and other operating 
parameters. 

Accelerations of point Pj^ in the Xq-Zq and X^-Yq planes are shown, respec- 
tively, in figures 4(e) and (f). These must be integrated over the volume of the nutator 
at a specific value of a?t usii^ equations (5) to obtain the moments about the inertial 
axes. These are then added vectorially to get the total inertial moment that must be re- 
acted by the bearing system that supports the nutator. 

Equations (2) to (5) can also be used to analyze squeeze film and distortion effects in 
hydrodynamic thrust bearings and face seals subjected to high-frequency nutations. 

Small -amplitude, high-frequency nutations develop in face seal runners. The extremely 
small film thicknesses in face seals coupled with the typically fragile structure of face 
seal runners lead one to believe that the inertial forces that result from nutating motions 
may significantly affect the performance of the seal. 


Inertial Moments 

Expressions for the inertial moments M , M , and M acting about the inertial 

0 ^0 ^0 

coordinate axes Xq, Yq, and Zq, respectively, are derived in appendix A for a nutating 
body of arbitrary shape. Equations (AlO) to (A12) are applicable at cot = 0. Equa- 
tions (A19) to (A21) are applicable at cot = t)’/ 2. Because the inertial moment must ob- 
viously be a rotating moment of constant magnitude for a constant frequency of nutation, 
its components can be calculated at any convenient value of cot and summed vectorially. 

The equations were derived for two values of cot as a check on their correctness. A 
comparison of equations (AlO) to (A12) with equations (A19) to (A21) shows their equiv- 
alence. In other words, for a circular disk of uniform thickness 
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M. 


= M. 


o;t=0 


'0 


M 


yo 


= M. 


o;t=0 


^0 


o;t=Tr/2 

a)t=7r/2 



a ±=0 



lwt=7r/2 


The derivation and calculation of inertial moments for a thick nutating plate are 
given in appendix B. As expected, M , the moment about the rotational axis (in the 

Xq-Zq plane) acting at the nutation frequency, is identically zero for a constant oj. At 
cot = 0 the inertial moment acts about the Zq axis, and at cot = ir/2 it acts about the 
Yq axis. Its rotational character is apparent from these relations. 

The magnitude of the inertial moment acting on the nutatirg plate used in the trans- 
mission design in reference 2 (the plate dimensions are given in appendix B) was cal- 
culated to be 7. Slxio"^ newton-centimeters (6. 66x10® lb-in. ). This exceeds by a factor 
of 28 the moment load imposed by the static axial reaction loads from the teeth of both 
the drive and reaction plates (ref. 2). It is quite obvious from these results that nutating 
plate drives are unusable at any operating conditions other than very low speeds. Even 
at very low speeds the magnitude of the reaction load and inertial moments on the nu- 
tating plate support bearing is likely to impose a severe bearing life problem. 

It is interesting to note from equation (B7) that, unless the plate thickness T is 
substantial relative to its radius r, the contribution of its thickness (as given by the 
first term of eq. (B7)) to the inertial moment is insignificant. For the plate dimensions 
used in the calculations, the first term in equation (B7) increased the inertial moment by 
approximately 0. 3 percent. Therefore, for relatively thin plates, the first term in equa- 
tion (B7) can be neglected. 


Comparison of Inertial Moments Due to Nutating and Simple Harmonic Motion 

Appendix C contains a derivation of the inertial moment developed in a body subject 
to harmonic oscillation. Equation (C7) gives the maximum inertial moment due to har- 
monic motion, which is developed in a circular plate. It can be compared with equa- 
tions (B4) and (B7). For a thin plate equation (C8) is compared with the second term in 
equations (B4) and (B7). Notice that the harmonic and nutating moment equations are 
identical except that equation (C8) contains the factor a (expressed in radians), while 
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the nutating body equations contain the factor tan o'(sin^o> + 1). 

A comparison of the magnitude of these terms as a function of a is given in table I, 
expressed as a percent error if equation (C7) is used, rather than the exact equation. 

The error is small at small values of a, but becomes appreciable as a increases. 
Thus it appears that the approximate inertial moment equation assuming simple har- 
monic oscillation can be used with good accuracy at small values of of. At larger values 
of a the exact nutation equations should be used. 

An additional value of the displacement, velocity, and acceleration equations 
(eqs. (2) to (4)) may lie in their use for determining local distortions in nutating bodies. 
Different points in a nutating body ejqierience varying accelerations at given times. This 
condition will result in varying stresses and local distortions. These may be significant 
in a machine element such as a face seal whose performance is a strong function of its 
geometry during operation. 


CONCLUSIONS 

The equations developed for the displacement, velocity, and acceleration of a point 
in a nutating body can be used for the design of nutating plate transmissions. The equa- 
tions developed for the calculation of the inertial moment of a nutating body of arbitrary 
shape, when applied to the plate of a previously designed nutating plate drive, indicate 
that this device is severely speed limited. The support bearing load due to the inertial 
moment was calculated to be 28 times as great as the load due to reaction forces. 

At small values of the nutation angle, the approximate formula for the inertial mo- 
ment, which is based on an assumption of simple harmonic oscillation, can be used with 
good accuracy. At nutation aisles of 10*^ and greater the error becomes appreciable, 
and the exact formula for nutating motion should be used. 

The equations for displacement, velocity, and acceleration should be useful in as- 
sessir^ squeeze film and local distortion effects in thrust bearings, and particularly face 
seals, in which nutating motions occur. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 8, 1973, 

501-24, 
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APPENDIX A 


DERIVATION OF INERTIAL MOMENTS FOR A BODY OF 


At o)t = 0 


ARBITRARY SHAPE AT oit = 0 AND cot = t/2 


“yo ‘*0^0 - Vo'P 

(Al) 

<Vo - Vo'P ■*'' 

(A2) 

■/" 

(A3) 

Xq = (cos Qf)xj - (sin o!)yj 

(A4) 

Yq = (cos Q!)yj + (sin a)x^ 

(A5) 

Zq = Zi 

(A6) 

_ co^ Sin a „ 

^0 2 

cos Of 

(A7) 

y. = Vsln„)x,.(“^^‘"^“)yj 
^ ^ \ COS a / 

(A8) 

- 2 2 
Zq = -(o) sin (3?)Zj 

(A9) 

•• ^ o)^ sin « „ ^ 

*0^0 = 2 

cos^ a 


2 2 2 3 

ZqXo = -(co sin a cos a)x^z^ + (co sin a)yj^Zj 
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^ 0^0 - Vo = " 


A 2 2 \ 

/I - sin or cos a\„ „ , . 

( ; ly^z^ + Csin 


Q! COS Cy)Xj^Zj 


cos a 
2 2 


/ 2 . /l - sin^a cos^«\ 

o> smo. ( ^ jyjz^ 

COS^Qt ' 


+ (sin a cos a)XjZ^ 


p dv 


2 = 


2 . 0/..2 „.„2 


2 sin^ a 2 




cos a 


yo^o = 


2 

2 sin a „2 , , 2 sin or 
y- + CO 

cos a 2 

cos Q! 


xjyi 


^0^0 " ^0^0 ~ tan 0!(1 + sin^of)yj + co^ tan^a(l - 2 cos^Qf)Xjyj + co^(sin a cos 

= - J [(sin a cos of)Xj + tan o!(sin^af + Dy^ + tan^o!(l - 2 cos^Q!)xjy jj ^ 


^0^0 “ sin^orcos 0 ')yjZj - (co^ sin^O')x^Zj 


y^Zg = -(co'^ sin a)x 


\ COS a / 


yi^i 


M 


= J" ^ sin a |(cos^Of)XjZj - [tan of(cos^Q! + 1) y^Zjj p dv 


At cot = 7 t/2 


Xq = Xj cos a + Zj sin or 


yo=yi 


Zq = Zj cos 0 ? - Xj sin a 


(AlO) 


a\x\] 

I dv 
(All) 

(A12) 

(A13) 

(A14) 

(A15) 
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2 . 

^ _ ;w sin a\^ 

O' 'i 1 

cos a 


(A16) 


Yq = -(w^ sin^Q!)yj 


(A17) 


t o 2 \ 
u) sin a\ 

1 

cos 0? / 


(A18) 


M. 


yo 


From equations (Al), (A13), (A15), (A16), and (A18) 

= - J' j^(sin a cos oi)xj + tan a(sin^<y + l)zj - tan^o!(l - 2 cos^O()Xj^Zjj p dv (A19) 


From equations (A2), (A13), (A14), (A16), and (A17) 


'•7 

t/y} 


M = / 0) sin Of 


2 2 

/ • \ il - sin Q! cos al „ _ 

(sin a cos aiXjYj - I J ^1^1 

cos^a 


p dv (A20) 


From equations (A3), (A14), (A15), (A17), and (A18) 


M. 


= J' (j?‘ sin a |(cos^a)xj^Zj + |tan or(cos^Q! + dv (A21) 


11 



APPENDIX B 


DERIVATION AND CALCULATION OF INERTIAL MOMENTS 
FOR A THICK NUTATING PLATE 


From figxire 5 


xi=xi 


Yj = r sin ^ 

Zj = r cos 4 
dv = r dr d^ dx^ 


From equation (A 10) and at wt = 0, 


M = sin a 
^0 


2 2 

1 - sin a cos a 
cos^a 


r2n rT/2 

Jq Jt. J-T /2 


r sin ^ cos ^ dxj^ dr d^ 


+ to 


1.2 r^o 2 

p sin a cos ® J J ^ 


Both integrals vanish so that 


M =0 
^0 


From equation (All) and at tot = 0 


(Bl) 


(B2) 


M = -io‘ 
^0 


! . /'2tt rv ^T/2 2 

p sin a cos a / / / x^r dx, dr 

^r. y-T/2 ^ ^ 


d| 


-co^p tan a (sin^a + 1) f f ^ 

'^0 Jt. Jj] 


27 t fT rT/2 

^ I r'' sin"^ dx, dr df 
j /-T/2 1 


-to 


! 2 2 /"r /'T/2 » 

p tan a(l - 2 cos J ^ J fY/2 ^ 


(B3) 
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M = -iroo^pT sin a cos a 
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2 2 
- TTO) pT tan o!(sin a + 1) 


..4 ^4‘ 


(B4) 


From equation (A 12) at 0)t = 0 


2 2 /*2jt /*T/2 2 

sin a cos a / / / x.r cos | dx- dr d^ 

^0 •'O Jr. At/2 ^ ^ 


9 2 /* /* 3 

tan a(cos'“a! + y y 2 ^ ^ ^ 


'0 *^r, ^-T/2 


(B5) 


Both integrals vanish, so that 


=0 

^0 


(B6) 


as expected if w is constant. From equations (A 19) to (A21), it can be shown that, 
at ^t ■ n/2 , 


M = -irat^pT sin a cos a 
^0 


12 


2 2 
JTO) pT tan Q!(sin a + 1) 


V4 A 

T - r . 
O 1 


(B7) 


M, =0 
^0 


=0 

^0 


(B8) 

(B9) 


As expected 


M, 


= M„ 


O)t=0 


^0 


wt=ix/2 


M„ =M 


lcut=0 nctJt=7r/2 


M„ 


= M. 


x; 


= 0 


la;t=0 


(i)t=v/2 
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At any value of wt the inertial moment Mj will be 

(BW) 

From these results it is clear that the inertial moment generated by the nutating plate 
acts about the axis that is in the YqZq plane and rotates with the nutation frequency w. 
For the specific plate used in the transmission design in reference 2; r^ = 45. 7 
centimeter (18 in.), r. = 20. 3 centimeter (8 in.), T = 4. 65 centimeter (1.83 in.), tt) r 587 
radians per second, o = 10°, and p= 7.9x10"^ N-secVcm^ (7.34x10”“^ lb- secVin.'^). 

With these values the inertial moment Mj becomes (from eqs. (B2), (B4), and (BIO)) 
7. 51x10 newton centimeters (6.66x10® lb* in.). 
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APPENDIX C 


INERTIAL MOMENT FOR A HARMONICALLY OSCILLATING PLATE 

From figure 5, if the plate executes an oscillation of amplitude 2 a at frequency co 
about the Zq or Z axis 


Tj = Of cos o)t 


77 = -aw sin (jjt 


rj = -aw cos 0)t 


The inertial moment due to harmonic motion will then be 


^h ~ ^ZZ^ 


The maximum moment occurs at a;t = 0 and tt where 


^ X 17 17 oico 

^max 




M ,,!!£T2-!r('r4 ,4^^ 

^max 4 L'® ^ ' 3'° 


where a is given in radians. 
For a thin plate 


M 


= ^ / 4 _ ^4\ 

^max 4 \ ® 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6) 

(C7) 


(C8) 


Again, for the specific plate used in the transmission design in reference 2 (dimen- 
sions and operating conditions given in appendix B), =7. 21x10*^ newton - 

centimeters (6. 39X10® lb- in. ). 
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